Nodal Structure of Superconductors with Time-Reversal Invariance 

and Z2 Topological Number 
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A topological argument is presented for nodal structures of superconducting states with time- 
reversal invariance. A generic Hamiltonian which describes a quasiparticle in superconducting states 
with time-reversal invariance is derived, and it is shown that only line nodes are topologically stable 
in single-band descriptions of superconductivity. Using the time-reversal symmetry, we introduce 
a real structure and define topological numbers of line nodes. Stability of line nodes is ensured by 
conservation of the topological numbers. Line nodes in high-Tc materials, the polar state in p-wave 
paring and mixed singlet-triplet superconducting states are examined in detail. 
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I. INTRODUCTION 

Nontrivial nodal structures are one of the most no- 
ticeable features of unconventional superconducting gap 
functions. Conventionally, the nodal structures have 
been investigated by power law behaviours of the tem- 
perature dependences of the specific heat, the NMR re- 
laxation rates and so on [1, 2]. Although the power law 
behaviours give a hint of the nodal structures, they can 
not give the details such as the number and the position 
of nodes. Recently, it has been increasingly clear that 
the angular-controlled measurements of thermal conduc- 
tivity and specific heat in a vortex state are useful probes 
to determine the details of the nodal structures [3-5]. 
Various unexpected nodal structures of superconductors 
have been found in this method [6-20]. This makes ur- 
gent the need for a better theoretical understanding of 
nodal structures of superconductors. 

The purpose of this paper is to reveal theoretically a 
generic nodal structure of superconducting states with 
time-reversal invariance. We use two different methods to 
investigate the nodal structures. The first one is based on 
a universal Hamiltonian we will construct, and the other 
is an argument using topological numbers in momentum 
space. Unless stated explicitly, our results do not rely on 
any specific symmetry of superconducting states except 
the time-reversal symmetry. 

Topological arguments in momentum space are power- 
ful tools to investigate a nonperturbative aspect of quan- 
tum theories [21-24]. For superconducting states with 
broken time-reversal invariance, a topological character- 
ization of the nodal structures was given in Refs. [25-27]. 
On the other hand, we will present here a topological 
characterization of nodal structures of superconducting 
states with time-reversal invariance. We will introduce 
novel topological numbers and discuss stability of the 
nodal structures. 



The organization of this paper is as follows. In Sec. II, 
we derive a generic Hamiltonian describing a quasipar- 
ticle in superconducting states with time-reversal invari- 
ance. The quasiparticle spectra are given and their nodal 
structures are discussed. In Sec. HI, we examine stabil- 
ity of line nodes for high-Tc materials, the polar state in 
p-wave paring, and mixed singlet-triplet states, respec- 
tively. In Sec. IV, a real structure is introduced and novel 
Z2 topological numbers of line nodes are defined. We cal- 
culate the Z2 topological numbers for high-Tc materials, 
the polar state in p-wave paring, and mixed singlet-triplet 
states. It is shown that the existence of the topological 
numbers ensures stability of the line nodes. Conclusions 
and discussions are given in Sec.V. 



II. HAMILTONIAN WITH TIME-REVERSAL 
SYMMETRY AND NODAL STRUCTURE 

We derive here a generic Hamiltonian which describes 
a quasiparticle in superconducting states with time- 
reversal invariance. Let us consider a single-band descrip- 
tion of superconducting state and start with the following 
Bogoliubov-de Gennes type Hamiltonian, 



i/-i^c(fe)tif(fc)c(fc). 



(2.1) 



where c(fc) is a four component notation of electron cre- 
ation and annihilation operators in momentum space. 



c(fc) 



c,(-fc)t 



ia = ±) 



(2.2) 



and H{k) is a 4 x 4 Hermitian matrix which will be de- 
termined below. Since c(fc) satisfies 



c(-fc)* = rc(fc), r 



1 

1 



(2.3) 



we can set the matrix H{k) to satisfy the following equa- 
tion. 
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TH{k)T = -H{-ky 



(2.4) 
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We also demand that the Hamihonian is time-reversal 
invariant. The time-reversal operation T is defined as 

Tc{k) = ec{-kr, e=(^y (2.5) 

where cxi's (i = 1,2,3) are the Pauli matrices. Time- 
reversal invariance of H implies 

eH{k)e-^ = H{-ky. (2.6) 

We solve Eqs.(2.4) and (2.6) to obtain a generic Hamil- 
tonian of a quasiparticle in superconducting states with 
timc-rcvcrsal invariance. 

To solve Eq.(2.6), we rewrite the 16 components of 
H{k) in terms of the identity matrix, 5 Dirac matrices 
r„ and 10 commutators r"*" = [ra,Tb]/{2i) [28], 

5 O 

H{k)=ho{k)l + J2ha{k)r,+ J2 hab{k)Tab.{2.7) 

a=l a<6=l 

Here ho{k), /ia(fe)'s and /iab(A;)'s are real functions of k. 
We choose the following representation of Ta, 

^1 = ( 1 ) ' ^^={1 -l)' 

/ -iai \ f -1(72 \ 

^'-[za^ j' ^'=[ia2 J' 

In this representation, we obtain 

er„e-^ = r:, QVabQ-' = -K^,. (2.9) 

Therefore, Eq.(2.6) is satisfied if we have 

ho{k) = ho{-k), ha{k) = ha{-k), (2.10) 

and 

habik) = -habi-k). (2.11) 

In addition to Eq.(2.6), we impose Eq.(2.4) on Eq.(2.7). 
Using r = Fi and the commutation relation of r^'s, we 
find that the following 8 functions should be identically 
zero, 

hoik) = hi{k) = h^{k) = /i5(fc) = 0, 

/ii3(fe) = /ii5(fe) = h2i{k) = h35{k) = 0. (2.12) 

Therefore, the Hamiltonian contains 16-8=8 independent 
functions. It is convenient to introduce a new notation 
of the remaining functions, 

h2{k) = e{k), hi{k) = -V'(fc), 
(/i45(fc), hii{k),hM{k)) = g{k), 
{h25{k), hi2{k),-h23{k)) = d{k). (2.13) 



All of them arc real functions. Here e(fc) and tp{k) are 
even functions of fc, and g(k) and d{k) are odd functions 
of fe. In terms of them, H{k) is written as 

where A(fe) is defined by 

l^{k) = ixp{k)(j2 + id{k) ■ cTa2. (2.15) 

Now physical meanings of these functions are evident. 
The function e(fe) is a band energy of electrons mea- 
sured relative to the chemical potential ji, and A(fc) is 
a gap function of a superconducting state. (■0(fe) and 
d{k) represent the spin-singlet and spin-triplet gaps, re- 
spectively). The fimction g{k) is a parity breaking term 
in the normal state. For example, the Rashba term is 
represented by this. 

The quasiparticle spectra E{k) in the superconducting 
state are given by the eigenvalues of H{k). The eigenval- 
ues can be obtained straightforwardly, and the resultant 
spectra E{k) are 

^(fc)= 

±yj e2 + + g2 + d2 ± 2^J{eg + i^df + {g x df 
= ±E±{k). (2.16) 

Zeros of E± (fe) determine the nodal structure of the su- 
perconducting state. First, it can be shown that E+{k) 
has no zero in general; For E^(k) to be zero, we have (at 
least) 

e{kf + iP{kf + g{kf + d{kf = 0. (2.17) 

This leads to 8 conditions, 

e(fe) = V(fe) = g{k) = d{k) = 0, (2.18) 

which can not be met generally in three dimensional mo- 
mentum space. Even if E+{k) has an accidental zero at 
some fe, we can easily remove this by a small deforma- 
tion of the Hamiltonian. Therefore, the zero of E^ (fe) is 
topologically unstable. 

Let us now consider the condition E-{k) = 0. The 
Hermitian property of H{k) ensures that the eigenvalue 
E- (fe) is real. Therefore, we have the inequality, 

e^+^^ + g^ + d''> 2V(e9 + V'd)' + (9Xd)2. (2.19) 

E^ (fe) is zero when the equality in (2.19) is attained. We 
rewritten this as 

+ V' + + d^ 

> 2 + -^;2rf2 g2fl2 g2^2 -(g.d- Cip)^ .{2.20) 

With fixed e^, t/j'^, g^ and d^, we can maximize the right 
hand side of (2.20) at g ■ d = etp. Then if we assume 
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g - d = 



ei/j, the equality in (2.20) is attained if we have 
= V'^ +9^- Therefore, E_{k) becomes zero when 



g{k) ■ d{k) ^ e{k)i,{k), 
e{kf + d{kf = yj{kf ^ 



9{k? 



(2.21) 
(2.22) 



We can show that these two equations are the neces- 
sary and sufficient condition for E^{k) ~ 0; If the equal- 
ity in (2.20) is attained under a condition different from 
Eqs.(2.21) and (2.22), the right hand side of Eq.(2.20) ex- 
ceeds the left hand side by imposing Eq.(2.21). This con- 
tradicts the inequality (2.20). Therefore, only Eqs.(2.21) 
and (2.22) lead to E_{k) = 0. 

The two equations (2.21) and (2.22) define two surfaces 
in three dimensional momentum space. They intersect in 
a line in general, and the intersection line gives zeros of 
E^{k). Small deformations of the Hamiltonian change 
the surfaces slightly, but the intersection line does not 
vanish. This means that for superconductors described by 
H{k), we have a line node in general and the line node is 
topologically stable. On the other hand, a point node is 
accidental and it can be removed by a small deformation 
of the Hamiltonian. 

It is interesting to compare this result with that of 
group theoretical analyses. Superconducting states are 
known to be classified by using a group theoretical 
technique and the generalized Ginzburg-Landau theory 
[2, 29-33]. The classification shows that point nodes ex- 
ist in a superconducting state with time-reversal invari- 
ance. For example. Table VI in Ref.[2] shows that the 
gap function in the representation of the tetragonal 
group has the form 



di (^k'j - — - X ky ykj; 



(2.23) 



This gap function preserves the time-reversal symmetry, 
and it has two point nodes on kz axis. This result does 
not contradict ours. Since the existence of the point 
nodes are due to a symmetry property of the F^ rep- 
resentation, we can remove these nodes by a small per- 
turbation which are not the F^ representation. Indeed, 
they are removed by the deformation 



d{k) ^ d{k) + rjzkz, {i] « 1) 



(2.24) 



where rizk^ is the Fj~ representation of D4h. Gener- 
ally, the group theoretical method does not ensure the 
stability of the nodal structure against any perturba- 
tions which change the representation of the group. On 
the other hand, our results are robust against any time- 
reversal invariant perturbations. 



III. EXAMPLES 

Using Eqs. (2.21) and (2.22), we can examine stability 
of line nodes. Here we consider three examples, 1) high- 
Tc materials, 2) the polar state in p-wave pairing, and 3) 
mixed singlet-triplet superconducting states. 



A. high-Tc materials 

In high-Tc materials, the parity is conserved, thus we 
have 



g{k) = d{k) = 0. 



(3.1) 



The band energy e{k) mainly depends on k^ and ky, and 
the Fermi surface given by e(fc) = is two dimensional. 
The gap function is 



m^^o{ki-ki). 



(3.2) 



In this state, E+{k) = E-{k), and we have four line 
nodes on the Fermi surface at fcp = (ifcp, ±fcF)/\/2- (See 
Fig.l.) We will show that the line nodes are stable. We 



\|/ = 




\|/ = 

(kY,kp)/{2 
k. 



e = 



FIG. 1: Line nodes in high-Tc superconductors. 

perturb the parameters of the Hamiltonian as 

+ Se, 



e - 
9 



5g, d 



E - 



d = 5d. 



(3.3) 



First, consider the case of Sg ^ 6d = 0. In this case, the 
perturbation does not break the parity symmetry. Now 
Eqs.(2.21) and (2.22) reduce to 



e{k) + de{k) ^ 0, ij{k) + Si;{k) ^ 0. 



(3.4) 



These equations can be satisfied by a small change of 
the position of the line nodes fc = fco ^ ^ = '^'o + ^^o- 
Indeed, expanding Eq.(3.4) around fe = fco, we have 



Ve(fco) • Sko - 
Vi^iko) ■ Sko 



Seiko) : 
- Sij{ko) 



0, 



0. 



(3.5) 



(Here we have used e(fco) — V'(^o) = 0.) For any Se 
and Sip, there exists Sko which satisfies these equations. 
Therefore, the perturbation moves the position of the line 
nodes, but it does not remove the line nodes. 

When Sg and Sd are not zero, Eqs. (2.21) and (2.22) 
become 

(e -I- Se) ■ {^P + S^) = Sg ■ Sd. (3.6) 
(e + Se)^ -i^b + S^jj)^ = Sg^ -Sd^, (3.7) 

These two equations give two curves in the {ip+Sip, e+Se) 
plane. As is shown in Fig. 2(a), these two curves always 
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have two intersection points near the origin. Therefore, 
we have two different conditions of Une nodes correspond- 
ing to the two intersection points. This impUes that each 
hne node sphts into two after the perturbation as illus- 
trated in Fig. 2(b). Although the number of line nodes 
becomes eight, the line nodes do not vanish by the per- 
turbation. 



£-|-5£ ky 




(a) (b) 

FIG. 2: (a) Two curves given by Eqs.(3.6) (the solid lines) and 
(3.7) (the dotted lines) . Here we show the case of 6g^ — 5d? > 
and Sg ■ 5d > 0. (b) A schematic figure of the line nodes in 
high-Tc superconductors with a parity breaking perturbation. 



These line nodes vanish if the time-reversal symmetry 
is broken. For example, by deforming V'(fc) 

ip{k) '0(fc) + iV: {v — const.), (3.8) 

we can completely remove the line nodes. 



B. polar state in p-wave paring 

Here we consider the polar state in p-wave paring. The 
polar state is given by 

e(fc) = — - M, d{k) = Aozk,, (3.9) 
Zm 

and 

^{k)=g{k)^0. (3.10) 

It is a solution of the Ginzburg-Landau theory neglecting 
the spin-orbit coupling [1]. Obviously, the gap has a line 
node on the equator. See Fig. 3. 

The line node is unstable, because, for example, the 
following deformation of the gap function removes the 
line node completely, 

d{k)-^ d{k)+7]{xk^ + yk,), (3.11) 

where 77 is a small real number. This line node is not 
given by an intersection line of Eqs.(2.21) and (2.22), 
thus it is accidental and topologically unstable. 



ky 



kz 

FIG. 3: A line node of the polar state in p-wave paring. 



C. mixed singlet-triplet states 

For noncentrosymmetric materials, the parity symme- 
try is broken in the normal states. In the presence of the 
spin-orbit interaction, the absence of inversion symme- 
try give rise to parity breaking terms g{k) in the normal 
states. In such systems, singlet and triplet parings can 
be mixed in the superconducting states [34] . We examine 
stability of line nodes in the mixed singlet-triplet super- 
conducting states. 

As a concrete example, consider a s -\-p superconduct- 
ing state proposed to explain the superconducting state 
in CePtaSi. In this material, the existence of line nodes 
was reported experimentally [12]. In CePtaSi, the ab- 
sence of inversion symmetry give rise to the Rashba in- 
teraction, 

aik) =a^{-xky+yk^), (3.12) 

where a is a real number. The following gap function was 
proposed to explain the line node in the superconducting 
state [35], 

ilj{k) = ^-0, d{k) = Aoi~xky + yk,). (3.13) 

Here V^o and Ap are real numbers. As is shown in 
Ref. [35], one can show that under a suitable choice of 
a, and Aq, E-{k) has two line nodes. {E+{k) has 
no node.) We show a schematic picture of the nodes 
in Fig. 4. As in the high-Tc case, it can be shown that 
the line nodes are stable against any small time-reversal 
invariant perturbations; The line nodes are intersection 
lines of the surfaces given by Eqs.(2.21) and (2.22) in 
momentum space. Under a small perturbation, 

e— >e + (5e, ip ^ + SiJ;, 

g^g + Sg, d-^d + 6d, (3.14) 

these surfaces moves slightly, but the intersection lines 
do not vanish unless the lines shrink to points. 
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where S is defined by 



FIG. 4: A schematic figure of line nodes on tiie Fermi surface 
of CePtaSi [35]. Tiie spiiere and tiie solid lines denote the 
Fermi surface and the line nodes, respectively. 



In a similar manner, we can show generally that 
line nodes are topologically stable in mixed singlet- 
triplet superconducting states with time-reversal invari- 
ance. For example, the line nodes proposed for Li2Pd3B 
and Li2Pt3B [36] are topologically stable. 



IV. Z2 TOPOLOGICAL NUMBER 

In the previous sections, we have seen that line nodes 
can be stable in superconducting states with time- 
reversal invariance. In this section, we show that the sta- 
bility is explained by the existence of topological numbers 
defined by wavefunctions of quasiparticles. The topolog- 
ical numbers are closely related to a real structure of the 
Hamiltonian H{k). Although H{k) is not a real matrix 
in general, it has another real structure due to the time- 
reversal symmetry. 



A. Real structure and Z2 topological number 

1. Real structure 

In order to find a real structure, it is convenient to 
consider the following 8x8 Hamiltonian, 



n{k) 



( H{k) 







\^ H{-k) 
and its eigenvalue equation, 

n{k)^{k) = E{k)^(k) 



(4.1) 



(4.2) 



The eigenvalue E{k) is the same as that of H{k). 

Using the time- reversal symmetry of H{k), we can 
show 



e 

-e 



= 1. 



(4.4) 



The relation (4.3) implies a real structure of 7^(fc). To 
see this, we introduce an operator JC = 'E.K where 
K is the ordinary complex-conjugation operator. The 
new operator obeys K,^ = 1, then we can think of it 
as a new complex-conjugation operator. The extended 
Hamiltonian 7i(fe) is "real" in terms of the new complex- 
conjugate operation. 



JCH{k)JC-^ = n{k). 



(4.5) 



Now we can impose the reality condition on the eigen- 
f unction 'J'(fc), 



^{k) =/C*(fc)(= S«'(fe)' 
Using an eigenfunction of H{k), 

H{k)u{k) = E{k)u{k), 



(4.6) 



(4.7) 



we find the following two independent eigenfunctions of 
Ti.{k) which satisfy the reality condition. 



and 



^i(fe) = -2 
Since u{k) is rewritten as 



u{k) 

-eu{ky 



u{k) 

eu{ky 



u{k) 




4'R(fe)-Hi*i(fc), 



(4.8) 



(4.9) 



(4.10) 



these two wavefunctions correspond to the "real" and 
"imaginary" parts of u{k), respectively. If u{k) is nor- 
malized as 



u{kyu{k) = 1, 
the normalization of \E'yi(fc) is given by 

^'A(fc)^*s(fc) - ^Sab, 

where A. B = R, I. 



(4.11) 



(4.12) 



2. Z2 topological number 

To define topological numbers of a line node, let us 
consider first an infinitesimal circle around the line 
node and solve the eigenequation H{k)u{k) = E{k)u{k) 
on S^. (See Fig. 5.) Because of the gauge freedom 



En{k)E-^ ^n{k)*, 



(4.3) 



u(k) 



u{k), 



(4.13) 
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u{k) is not determined uniquely. We have to fix the 
gauge. Generally, a gauge fixed solution has a singularity 
on S^, thus it can not be well-defined on the entire S^. 
Two or more solutions with different gauge fixings are 
needed to cover the entire S^. We consider a couple of 
solutions u^^)(A;) and u'^^^{k), and demand that the first 
component of u'^^^{k) and the second one of u^'^'>{k) are 
real. Because of these different demands, the solutions 



(fc) and u^^-* (fc) have different singularities from each 



other. Therefore, we can cover the entire 5*^ by using 
them. 

Let us next consider the Kramers doublet of u'^^'>{k), 
eu(2)(-fc)*. By definition, the Kr amers doublet has the 
real first component as well as u^^'^^k). Therefore the 
following two possibility arises, 

1. eii(2)(_fc)* = (or -u^^Hk)), 

2. ew(2)(_fc)* ^±M(i)(fe). 

Here we concentrate on the latter case, since only this is 
relevant to the topological numbers of the line node [47] . 
In this case, 6m^^)(— A;)* has a different singularity from 
M^^)(fc), thus we can use Qu^^\—k)* instead of u'^-'(fc) 
to cover the entire region of S^. After all, the following 
two solutions with the real first component are obtained, 

uW(fe) = M(i)(fc), w(-)(fc) = eu(2)(-fe)*. (4.14) 

If the solutions u'^^^k) and u^^\k) are nonsingular on 
U+ and U-, respectively, is given by the union of C/+ 
and [/_, 



(4.15) 




= A 



II 



u. 



FIG. 5: around a line node. is covered by U+ and U- 
The overlap of U+ and U- consists of the regions I and II. 



Now it is convenient to introduce the eigenfunctions 
of Ti{k) instead of M^'^-'(fc) (a — ±). The eigenfunctions 
of 7i(fe) with the reality condition (4.6) are constructed 



from M^""^ (fc) as 



*{"'(fc) = - 



u(^)(fc) 

-eu'-'^^k)* 

i_ f u('^)(fe) 
2 eu('^)(fc)' 



(4.16) 



It is easily find that the eigenfunctions 4'|^''(fc)'s {A 
R, I) are regular in [/o-. 



On the overlap [/+nf7_ 



both ^'^^(A;) and ' (k) are 



regular and they are related by the transition function 

TAB{k), 



f^+)(fe)=T^B(fc)*r(fe) 



(4.17) 



The reality condition (4.6) and the normalization (4.12) 
imply that the transition function T(fc)(= TAsik)) is an 
element of 0(2). Moreover, T{k) reduces to 



T{k) = 



_ ( tn{k) 



ti{k) 



iR,(fe) =ti(fc) =±1,(4.18) 



since ^'p'^''(fc)'s have the real first components but 

¥^''\kys have the imaginary ones. Therefore, ^'^+'(fc) 

(or 'i'['^\k)) is identified with ^(^"^(fc) (*{~^(fe)) on 
C/+ n U-. As we will show immediately, the transition 
function determines the global topology of the wavefunc- 
tions. 

For the sake of simplicity, we assume in the following 
that both U+ and U- are connected regions as illustrated 
in Fig. 5. In this case, the overlap [/+ n C/- consists of 
two regions I and II. The generalization to other cases in 
which C/+ and U- consist of many disconnected segments 
is straightforward. 

Let us examine the topology of ^'^'(fc). On S\ it 
has two transition functions tR(fci) and tR(fcii) where 
fci and kii are arbitrary points on the regions I and II, 
respectively. If t^{ki) = tR(fcn), ^'|^^(fe) are glued to 
^t^l^ \k) without twisting, thus we have a trivial topol- 
ogy similar to a cylinder. (See Fig 6.) This configura- 
tion is a contractable loop in the Hilbert space, there- 
fore nothing interesting happens. On the other hand, if 

tR(fci) = -iR(fcii), *R'^''(fc) and *R~'(fc) are glued to- 
gether with twisting, thus we have a nontrivial topology 
similar to the Mobius strip [37] . Since this configuration 
is a noncontractable loop around the line node, the line 
node can be considered as a kind of topological defect 
(vortex). Therefore, in this case the line node is sta- 
ble against a continuous deformation of the Hamiltonian. 
If we twist the wavefunction in twice, we have a trivial 
configuration again. Thus, the corresponding homotopy 
group is TTi = 

By this argument, we naturally introduce the following 
Z2 topological number /r. 



In (iR(fei)tR(feii)) 



(mod.2) 



It takes 



0, for iR(fci) = tR(fcii), 

1, for tn{ki) = -tn{ku). 



(4.19) 



(4.20) 



The line node is topologically stable when /r = 1 
(mod.2), and the stability is explained by the conser- 
vation of the topological number. 
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II 



(+) 



(-) 





II 



FIG. 6: Possible configurations around a line node, (a) A 
contractable loop, (b) A noncontractable loop. 



In a similar manner, another Z2 topological number 
can be defined as, 



h = 



In (ti(fci)ii(fcii)) 



(mod. 2). 



(4.21) 



Since ti{k) is equal to iR(fc), it is easily seen that /j = /r. 

Unless the time-reversal symmetry is broken, the two 
topological numbers /r and /i are conserved indepen- 
dently. However, once the time-reversal invariance is 
lost, they are not conserved. The reality of the Hamil- 
tonian (4.6) is lost and 5'R(fc) cannot be distinguished 
from ^'i(A;). Although the summation of /r and /j may 
be conserved, it is identically zero 



Ik + /i = (mod.2) 



(4.22) 



Therefore, the line node loses the topological stability 
and it disappears in general. 

It can be shown that the topological numbers are 
gauge invariant if the gauge transformation is nonsin- 
gular. Consider the following gauge transformation 



u{k) u'{k) = e*''Wii(fc) 



(4.23) 



where 9{k) is a nonsingular function of k. When 9(k) 
is nonsingular, the regions U+ and U- in which u^^'^'{k) 
and -*'(&) are nonsingular remain the same as before. 
The U{1) gauge transformation (4.23) induces the follow- 
ing S0(2) gauge transformation on \E'^^(fe), 



■^^^^'{k) = OAB{k)-^';'{k), 



where 



Oik) = OAB{k) 



cos 9{k) —sm6{k) 
sm6{k) cosd{k) 



(4.24) 



(4.25) 



The new transformation function T'{k) between 



\I'|^'''(fc)'s is given by 



T'{k) = 0{k)T{k)0{ky 



(4.26) 



Because the 0(2) matrices Oik) and T{k) commute with 
each other, we obtain 



r'(fc) ^ T{k). 



(4.27) 



Therefore, the topological numbers /r and /j defined by 
r(A;) are gauge invariant. 



B. high-Tc materials 

The Hamiltonian of high-Tc superconductors is 
e(fc) iip{k)a2 



H{k)^ 



-iip{k)a2 -e(fc) 



(4.28) 



This superconducting state is unitary, thus the quasipar- 
ticle energies E±{k) are degenerate in the entire momen- 
tum space, 

E±{k) = E{k) = y/e{k)^ +tp{k)^. (4.29) 

Since the parity is conserved and H{k) is real, H{k) com- 
mutates with 8, 



eH{k) ^ H{k)e. 



(4.30) 



Therefore, the eigenfunctions Ui{k) [i = 1,2) with E{k) 
can be the eigenfunction of O simultaneously, 

eMi(fc) = mi(fc), Qu2{k.) ^ -iu2{k). (4.31) 

Consider a line node shown in Fig. 7. The eigenfunctions 

\(/ = 
\[/<0 \|/>0 




8>0 



FIG. 7: around a line node in high-Tc superconductors 



regular in C/+ are given by 

u']^\k) 



2y/E{E'- 



ui+\k)- 



2y/E{E + e) 



/ {E + e)\ 

i{E + e) 
—it] J 
\ ^ J 

( E + e 
-i{E + e) 
iip 



, (4.32) 



and those regular in [/_ are 




(4.33) 



In a similar manner as tA^k) and I a {-A — R, I) in the 
previous subsection, the transition functions taAi^) and 
the topological number la a (a = 1, 2, A = R, I) are con- 
structed from Ua\k,) (a = 1,2, cr = ±). The transition 
functions become 



ti^{k) = tn{k) = sgmjj{k), 
hnik) = t2i{k) = sgml>{k). 



(4.34) 



Because the transition functions have different signatures 
between regions I and II in Fig. 7, the topological numbers 
calculated on 5*^ are given by 



'IR 



1, 



i2R 



21 



1. 



(4.35) 



As long as the parity is conserved, the eigenfunctions of 
H{k) stay the eigenfunctions of Q. Therefore, Iia and 
I2A {A ~ R, I) are not mixed, and the stability of the 
line node is ensured by conservation of these topological 
numbers. 

If the parity is broken, the degeneracy between E^{k) 
and E-{k) is resolved and only E-{k) has line nodes. 
When the perturbation is small enough, the topologi- 
cal numbers /r and Ii defined by the eigenfunction with 
E^{k) are calculated on the same as 



IR 



2R 



0, 



21 



0. 



(4.36) 



This is consistent with the splitting of the line nodes 
shown in Sec. Ill A; The line node in Fig. 7 is divided into 
two as is shown in Sec. Ill A, and each divided line node 
has the topological numbers /r = /i = 1. Because en- 
closes both the divided line nodes, the topological num- 
bers calculated on become /r = /i = 1 + 1 = which 
is the same as Eq.(4.36). 



C. polar state in p-wave paring 

The Hamiltonian of the polar state in p-wave paring is 

= ( Jtl ) (4-37) 



The quasiparticle energies E±{k) = E{k) are given by 
E{k) = v/e(fc)2 + d3(fe)2. (4.38) 




FIG. 8: around the line node in the polar state. 



Consider a circle around the line node in the equator. 
(See Fig. 8.) The circle is covered by the following 
solutions, 



4+'(fc) 



u'^-\k) 



4+\k) 



u[-\k) = 



1 


(Y) 


2y/E{E^ 
















1 






2^E{E- 




E-e 
\E-eJ 

( E + e 




(4.39) 



Using these solutions, we obtain the transition functions 
as follows, 

iiR(fc) = tii{k) = sgnd3(A;), 

t2R(fc) = t2i{k) = -sgnd3(fc). (4.40) 

From this, the topological numbers are calculated as 

/iR - hi = 1, /iR = hi = 1. (4.41) 

At a first glass, it seems to suggest that the line node has 
nontrivial topological numbers. However, on the con- 
trary to the high-Tc case, there is no quantum number 

which distinguishes u^*^' (fc) from Mj'^'' (fc) (cr = ±). There- 
fore, only the summations Iia + I2A — (^ = R, I) are 
preserved, and the line node has no nontrivial topolog- 
ical numbers. This result is consistent with the result 
in Sec. Ill B; The line node is unstable against a small 
perturbation of H{k). 

D. mixed singlet-triplet state 



The Hamiltonian of mixed singlet-triplet states is 



H{k) 



e{k)+g{k)-CT 
A(fc)t 



A(fc) 
-e{k) + g{k) ■ a* 



,(4.42) 
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where 

A(A;) = itp{k)a2 + id{k) ■ aa-i- (4.43) 

Here V(^) and d(fe) are real functions. As was shown in 
Ref.[38], d(fc) is proportional to g(fe) in mixed singlet- 
triplet states, 

d(fe) = cg(fc). (4.44) 

For the sake of simplicity, we assume c > in the follow- 
ing. -E-(fc) is given by 



= V(e(fe) - sgnV(fe)|g(fe)|)2 + (IV(fe)l - |d(fc)|)2. 

(4.45) 

The eigenfunction has the following two forms. 



u(-)(fc) = 





+ e-sgn-0|g|) 


( iE-+e 


- sgni^\g\) (j) \ 


\ -isgnip ( 


V|-|d|)a2<^; 




1 



y/2E_{E_ - e + sgmp\g\) 



X ( ^^gn^(l^|-|f)f ,y(4.46) 

where 0(fc) satisfies 

d{k) ■ cr(j){k) = -sgnV'(fe)|d(fe)|(/)(fe). (4.47) 
The first components of u^^^(A;) are real if we choose 

m = 

' i '^'"'^^ V (forsgnV^>0) 

-, ( i"^' , (for sgn^ < 0) ' 

(4.48) 

The transition functions obtained from 'u(+)(fe) and 
u(~^(fe) are 

tR(fe) = ti(fe) = sgnV'(fe) • sgn(|V(fe)| - |d(fe)|). (4.49) 

From these transition functions, we can calculate the 
topological numbers /r and 1\ of line nodes in mixed 
singlet-triplet superconducting states. For example, the 
straightforward calculation shows that both the line 
nodes of CeSisPt in Fig. 4 have nontrivial topological 
numbers, /r = 7i = 1. We can also show that the line 
nodes proposed for Li2Pd3B and Li2Pt3B [36] have non- 
trivial topological numbers. These results are consistent 
with the stability of the line nodes examined in Sec. Ill C. 

Note that the time-reversal symmetry is essential for 
the stability of the line nodes. For example, when '^{k) 



becomes complex, the time-reversal symmetry is broken. 
In this case, the quasiparticle spectra are given by 

E{k) 

= ± V(e(fe) ± l9(fe)l)' + (ReV'(fe) ± \d{k)\f + (ImV(fc))^ 

(4.50) 

If Im'!/;(fc) is a nonzero constant, the line nodes in the 
mixed singlet-triplet states vanish completely. 

V. CONCLUSIONS AND DISCUSSIONS 

(1) First we would like to outline the results of this 
paper. We examined topological stability of line nodes 
in superconducting states with time-reversal invariance. 
A generic Hamiltonian of time-reversal invariant super- 
conducting states was presented. We found that only 
line nodes are topologically stable in single-band descrip- 
tions of superconductivity. It was shown that line nodes 
in high-Tc materials and mixed singlet-triplet supercon- 
ducting states are topologically stable, while one in the 
polar state is not. Using the time-reversal symmetry, we 
introduced a real structure and defined topological 
numbers. Stability of line nodes was explained by con- 
servation of the topological numbers. 

(2) Besides the superconducting states examined in 
this paper, several superconductors such as CeCoIns [42], 
Celring [43] and Sr2Ru04 [7, 15, 16, 44] are believed 
to host line nodes. Among them, the gap function of 
CeCoIn5 is a rf-wave paring, thus its line nodes are topo- 
logically stable in a similar manner as high-Tc materi- 
als. On the other hand, the line nodes of Sr2Ru04 are 
not topologically stable since its superconducting state 
breaks the time-reversal symmetry [45] . 

(3) An analysis based on the KR theory [39] shows that 
point nodes have a Zi topological number in a nonrel- 
ativistic Fermi system with the time-reversal symmetry. 
However, our analysis in Sec. II indicates that point nodes 
have only the trivial topological number if the supercon- 
ducting state is described by single-band electrons. 

(4) If the superconductivity occurs in multi-bands of 
electrons, a topologically stable point node is possible to 
exist. In other words, if the time-reversal symmetry is 
not broken, the existence of a topological stable point 
node implies a multi-band superconductivity. 

(5) For mixed singlet-triplet superconducting states, 
the existence of the line nodes is not explained by the 
group theoretical method in Refs. [29-33]. The topologi- 
cal stability studied here is essential for the existence of 
these line nodes. 

(6) In addition to the time-reversal symmetry, if the 
party is conserved, the Hamiltonian becomes 

This is a real symmetric matrix, thus we have a natural 
real structure. Using this real structure, we can intro- 
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duce a Z2 topological number which is different from 
that given in Sec.IV as follows [40]. Let us consider the 
eigenfunction of H{k)u{k) = E{k)u{k), 

/ a(fc) \ 

uik) = ^[^j , (5.2) 
V d{k) J 

Here we can demand a{k), b{k), c{k) and b{k) to be 
real since H{k) is a real matrix. If we impose the nor- 
malization condition on ^(A;), we have a{k)'^ + biltf + 
c(fc)^ + d(fe)^ = 1, and by using the gauge freedom of 
the cigcnequation, u{k) can be identified with —u{k). 
Therefore, u{k) is given by an element of / [48]. 

Now consider a map determined by u{k) from an in- 
finitesimal around a line node into S'^ j Z2. The map 
is classified by t^xiS^ j Z^). Using the homotopy theory 
[41], we obtain 

7ri(5V-^2)=^o(-^2)-2^2. (5.3) 

Therefore, we have a .Z2 topological number correspond- 
ing to 'KxiS^ I Z'l). The line node is topologically stable if 
the map corresponds to the nontrivial element of ^2- A 
straightforward calculation shows that the map around 
a line node in high-Tc materials gives the nontrivial ele- 
ment. This is another explanation of topological stability 
of the line nodes in high-Tc materials. 

Although the construction of this topological number 



is easier than that given in Sec.IV, it is available only 
when the parity is conserved. 

(7) Wo have noticed that when a line node splits in two 
as was shown in Sec.IIIA, there remains a degenerate 
point between E+{k) and E+{k) near the fine node. The 
degenerate point is topologically stable since it is given by 
an intersection point of the three equations e{k)g{k) + 
ip{k)d{k) = in three dimensional momentum space. 
The degenerate point has a topological number similar to 
that of the chiral fermion in Ref.[21]. In order to describe 
the splitting of the line node in terms of the topological 
numbers given in this paper, we need to take into account 
the conservation law of this also. 

Note added.- After completing this work, the author 
noticed a preprint [46] by G. E. Volovik in which the 
topological stability of line nodes in high-Tj. superconduc- 
tors was discussed very recently. There is some overlap 
between his paper and Sec.IV B in this paper. His argu- 
ment was restricted to the case where the Hamiltonian is 
real, while the argument described here is not restricted. 
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